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ABSTRACT 

If gamma ray bursts are highly collimated, radiating into only a small fraction 
of the sky, the energy requirements of each event may be reduced by several 
(~ 4-6) orders of magnitude, and the event rate increased correspondingly. 
Extreme conditions in gamma ray bursters lead to highly relativistic motions 
(bulk Lorentz factors F ^ 100). This results in strong forward beaming of the 
emitted radiation in the observer's rest frame. Thus, all information on gamma 
ray bursts comes from those ejecta emitted in a narrow cone (opening angle 
~ 1 /r) pointing towards the observer. We are at present ignorant of whether 
there are ejecta outside that cone or not. 

The recent detection of longer wavelength transients following gamma ray 
bursts allows an empirical test of whether gamma ray bursts are well-coUimated 
jets or spherical fireballs. The bulk Lorentz factor of the burst ejecta will 
decrease with time after the event, as the ejecta sweep up the surrounding 
medium. Thus, radiation from the ejecta is beamed into an ever increasing solid 
angle as the burst remnant evolves. It follows that if gamma ray bursts are 
highly collimated, many more optical and radio transients should be observed 
without associated gamma rays than with them. Published supernova searches 
may contain enough data to test the most extreme models of gamma ray 
beaming. We close with a brief discussion of other possible consequences of 
beaming, including its effect on the evolution of burst remnants. 

Subject headings: Gamma rays: bursts — X-rays: general — 
ultraviolet: general — radio: general — ISM: jets and outflows 
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1. Introduction 

Relativistic expansion is a generic feature of fireball models for gamma ray bursts. This 
can be seen in several ways. First, the bursts are luminous (peak ~ in gamma rays 

for cosmological distances) and vary on millisecond time scales. The variability requires 
them to be small (size ^ a few x 10^^(1/100)2 cm [Woods & Loeb 1995]). This in turn 
guarantees that they exceed the Eddington limit on luminosity of an object in equilibrium, 
because the Eddington mass for L ~ lO^^L© comfortably exceeds the mass required to form 
a black hole in a region ^ 10^^ cm across. 

The observed gamma ray fluences imply total burst energies of order lO^^erg in the 
gamma ray range (20-2000keV) (e.g., Paczyhski & Rhoads 1993), assuming spherical 
symmetry. Together with the size limit, this implies energy densities ;^ lO^^erg/ cm^. 
This limit is sufficient to imply a large optical depth to electron-positron pair creation 
(77 ~^ e~^e~). Thus, the energy will not escape directly as electromagnetic radiation, but 
will be converted into a relativistic wind of pairs and baryons. The observed radiation is 
presumably generated later by interactions within such a wind (e.g., Paczynski & Xu 1994; 
Rees & Meszaros 1994) or with an ambient medium (Rees & Meszaros 1992; Meszaros & 
Rees 1993). The nonthermal gamma ray spectra provide additional support for this scenario 
(Goodman 1986; Paczyhski 1986; Fenimore, Epstein, & Ho 1993; Woods & Loeb 1995). 

If we are to understand the physical origin of the gamma ray bursts, we need to know 
about their distance, their luminosity, their frequency, and their environments. The recent 
detection of a variable optical counterpart to gamma ray burst 970508 (Bond 1997) and 
the detection of interstellar absorption fines at redshift z — 0.835 in its spectrum (Metzger 
et al 1997) clearly shows that at least some bursts are at cosmological distances, while 
the observed isotropy of the burst distribution on the sky (Briggs et al 1996) supports a 
cosmological origin for the overwhelming majority of the population. With the distance 
scale established, the largest remaining uncertainty in the burst energy is the solid angle 
into which the bursts radiate. This is also the dominant uncertainty in the burst event rate. 
If bursters beam their gamma rays into solid angle fl^, the burst energy scales as fl^/ (47r) 
and the event rate as Att/^Ij relative to the case of isotropic emission. 

The maximum plausible degree of collimation for electromagnetic radiation from a 
burster is opening angle (.y ~ l/F, where F is the bulk Lorentz factor of the radiating 
matter. This is the opening angle into which photons emitted isotropically in the rest 
frame of the radiating matter travel in the observer's frame (cf. Rybicki & Lightman 1979; 
Pacholczyk 1970). It is of course possible for the burst to be less colfimated, if the ejecta 
move into a cone of opening angle Cm > C7) but we still see only those photons emitted by 
matter in the smaller cone of opening angle ^7- 
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CoUimated jets are remarkably common in astrophysical sources. We observe them 
at small scales (protostars) and large scales (radio galaxies). The most widely accepted 
taxonomy of active galactic nuclei relies on orientation effects in accretion disk plus jet 
models to explain a variety of spectral features. Some sources have been observed with 
relativistic bipolar outflows (e.g., Galactic micro-quasars [Mirabel and Rodriguez 1994]), 
and these can show a marked asymmetry in apparent brightness between the approaching 
and receding jets that is well modelled as an effect of relativistic beaming (Mirabel & 
Rodriguez 1994; Bodo & Ghisellini 1995). Some possible mechanisms for gamma ray bursts 
allow naturally for beaming, through either very strong magnetic fields, accretion disks, or 
a combination of the two (e.g., the Blandford & Znajek [1977] mechanism). 

Lower energy transient sources are expected to follow gamma ray bursts as the fireball 
ejecta plow into the surrounding medium (Paczyhski & Rhoads 1993; Katz 1994; Meszaros 
& Rees 1997a; Waxman 1997; Wijers, Rees, & Meszaros 1997). The spectrum at a fixed 
time is generally modelled as a broken power law, resulting from synchrotron emission from 
a power law distribution of electron energies in a magnetic field at or reasonably near the 
equipartition value. The break in the spectrum shifts to lower frequencies as the burst 
remnant ages, primarily because the bulk Lorentz factor of the ejecta decreases, reducing the 
relativistic blue shift of the emitted spectrum. Generally, the break frequency is expected 
to decrease as a power law in time since the burst. While our imperfect understanding 
of relativistic, magnetized shocks leaves large uncertainties in these models, they are now 
observationally justified by the observed X-ray, optical, and radio counterparts to bursts 
970508 (lAU circulars 6654-6663) and 970228 (Wijers et al 1997, and references therein). 

Because the shift to lower frequencies accompanies the shift to lower bulk Lorentz 
factors, the minimum solid angle into which the transient can radiate increases with time. 
This leads directly to our proposed test for isotropy of gamma ray burst emission. If bursts 
are highly coUimated, the gamma rays will radiate into a small solid angle; the optical 
transient into a larger one; and the radio transient into a larger one still. Thus, we expect 
to see more optical transients than gamma ray bursts, and still more radio transients. 
On the other hand, if gamma ray bursts emit isotropically, we do not expect there to be 
optical transients unaccompanied by gamma ray bursts. The ratio of event rates for burst 
transients at two frequencies thus gives the ratio of the mean solid angle into which the 
burst transients radiate at those frequencies. We know the gamma ray burst rate well 
already, and within a few months should have a reasonable statistical sample of optical 
counterparts to observed gamma ray bursts. Establishing the total event rate for all optical 
transients with the characteristics of observed burst counterparts (whether gamma rays are 
seen or not) is a large but quite feasible task with present instruments. 
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2. Model-Independent Limits on Beaming 



The simplest model independent form of our test states simply that the ratio of event 
rates at two observed frequencies z/® i and z/0 2 should be N1/N2 = where Qi is the 

solid angle into which the flux is beamed at frequency and where iVj is the event rate 
at integrated over all fluxes. 

Unfortunately, we do not have the luxury of infinitely sensitive instruments. We 
therefore need to replace iVj with iVj, the rate of events at u^^i exceeding a flux detection 
threshold /min,i- Observed differences in Ni, N2 can then be explained either by different 
degrees of coUimation at different frequencies, or by insufficient sensitivity to detect some 
transients at one or another frequency. (The situation becomes more complicated if we are 
unable to select a sample of transients caused by a common physical mechanism.) 

To account for flux thresholds, consider the joint probability distribution for a 

burst to have angle-averaged fluxes /i, /2 at our two frequencies. ( "Angle- averaged" means 
that / = /Q/(47r) = L/ (AncP), where / is observed flux, L the source luminosity, and d the 
luminosity distance to the source.) Then our observations at frequency ^0^1 will detect a 
fraction Fi of all transients, where 



i^l = T^X / / p{h,f2)df2dh . 

47r ^/mi„,ini/47r Jo 

A similar equation gives F2, while the fraction of events seen at both frequencies is 

min [fli, 0,2) 



F\2 — 



47r 



X 



/ / P{hj2)df2dh 

•'/min,lf2l/47r -'/min,2fi2/47r 



(2) 



We then examine the ratio F12/F1: 



F12 
Fi 
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X 



min (fli, Q2) 



(3) 



The term in square brackets is < 1 because p is strictly non-negative. Also, 
min f22)/f2i < ^12/^1- Thus, ^12/^^1 > F12/F1 regardless of the details of the 
joint flux distribution p. Likewise, fli/fl2 > F12/F2. Although the Fj are deflned as 
fractions of the (unknown) total transient population, the ratios we care about can be 
expressed in terms of measurable event rates: F^/Fj — N^/Nj. Combining these results, we 
find 

(4) 



N12 Qi Ni 
— < — < — 

N2 ~ n2~ N12 



It is of course necessary that the flux thresholds /niin,i, /min,2 used to measure A^i2 be the 
same ones used to measure A^i and A^2- 
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If we can constrain p well, we may go beyond this analysis and actually estimate 
Qi/Q2- The resulting estimate will be sensitive to errors in p, however, and such errors 
will remain substantial until the lower frequency counterparts of gamma ray bursts are 
better studied. For now, we prefer to emphasize the inequalities in equation ^ which are 
model- independent. 

3. Expectations in an Illustrative Model 

We now turn our attention to modelling the expected transient event rate as a function 
of frequency. We are most interested in Tp^u^), the bulk Lorentz factor of the ejecta when 
the flux density peaks at observed frequency z/®. As with most predictions of relativistic 
fireball models, Tpi^u^) turns out to be a power law Tp cx over a large range of i/© and 
Tp. Pubhshed values of fi cover a substantial range, from /i = 1/4 (for the "impulsive 
fireball" model in Meszaros & Rees 1997a) to = 9/16 (from Paczyhski & Rhoads 1993). 
More secure quantitative predictions may be expected later, when uncertainties in the input 
physics of the models (particularly the electron energy spectrum) have been removed by 
confrontation with newly available counterpart observations (cf. Waxman 1997). 

The detection of radio emission from gamma ray burst 970508 (Frail et al 1997) 
at z > 0.835 suggests that rp(8.46GHz) ^ a few. Extrapolating boldly to optical and 
gamma-ray wavelengths, and using the relatively conservative scaling exponent jU = 1/4, we 
infer rp(R band) ^ 30, and Tpi^^y) = ^ 500. This Fq is consistent with the observational 
constraints of Woods & Loeb (1995). The maximum beaming for the gamma ray regime 
(taking F = 500) is then into 10^^ x An steradians. The optical regime (F = 30) gives 
2.8 X 10"'^ X Att steradians, while the radio (F = 2) gives 0.06 x 47r steradians (though 
the ( (X 1/r scaling is not very accurate at such low F). The precise values of F are less 
important than their ratios. 

The observed gamma ray burst rate from BATSE (the Burst and Transient Source 
Experiment aboard the Compton Gamma Ray Observatory satellite) is about 1 burst per 
day (Meegan et al 1996). The four bursts well localized by the BeppoSAX satellite (970111, 
970228, 970402, & 970508) have resulted in two probable optical counterpart detections. 
Thus, the rate of optical counterparts to observable classical gamma ray bursts is of order 
200 per year, or 1 per square degree per 200 years. We infer that if gamma ray bursts are 
maximally beamed, there should instead be of order 1 optical transient per square degree 
per year, and of order 1 radio transient per square degree per day. 

Whether the radio transients predicted under this extreme beaming scenario would be 
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visible is doubtful: The ratio of their peak radio flux to gamma ray fluence will be reduced 
by a factor [rp(radio)/rp(7)]^ ~ 10~^ relative to the isotropic burst case, assuming that all 
burst ejecta have about the same initial Lorentz factor. The reduction of flux would be 
less dramatic under scenarios where ejecta have a wide range of Lorentz factors, like the 
"layered jet" model (Meszaros & Rees 1997b; Wijers et al 1997) and "hypernova" model 
(Paczyhski 1997). In these cases, material ejected from the burster at comparatively low 
r might contribute substantially to afterglow flux at low energies and neghgibly at high 
energies. 

4. Observational Prospects 

The rate of gamma ray bursts is already fairly well known. We here consider the 
requirements for determining the rates of transient optical and radio events like that 
observed following gamma ray burst 970508. 

The peak optical emission of the 970508 transient was about R = 19.78 ± 0.05 
magnitude (37 ± 2 fiJy) (Mignoli et al 1997). The counterpart to burst 970228 was a little 
fainter (pa 17/iJy at V and 1 bands [Groot et al 1997]). The events lasted a few days each. 
We will assume that transients of this nature would be detected by daily observations to 
limiting 5(7 sensitivity i? 22. Such observations take a few minutes per fleld with a 1 
meter class telescope. The best fleld of view presently available at Kitt Peak National 
Observatory (1°, using the Mosaic CCD camera on the 0.9 meter telescope) would allow 
one to survey roughly 3 flelds per hour, or of order 20 square degrees per night, allowing for 
overheads. Thus, event rates of the order discussed in the extreme beaming scenario could 
be tested in ~ 10-20 nights on the telescope. 

Existing optical data may already be sufficient to apply a crude version of our test. 
Deep supernova searches by two groups (the High Redshift Supernova Search and the 
Supernova Cosmology Project) are ongoing and have now detected at least 37 (Schmidt et 
al 1996) and 28 (Dcustua et al 1996; Kim et al 1997) supernovae respectively. Their search 
strategy often includes observations separated by sufficiently short times (a few nights) that 
optical transients like the 970508 afterglow would be seen on multiple nights and confirmed 
as sources. We can estimate their total coverage as follows: The supernova rate is about 
34 per square degree per year for 21.3 < it! < 22.3 (Pain et al. 1996). The detected 65 
supernovae therefore translate to a coverage of about 2 degree^ year for supernovae. The 
effective coverage for gamma ray burst counterparts may be somewhat lower owing to the 
shorter time scale and hence reduced detection efficiency for these events. We will take the 
surveys to cover ~ 1 degree^ year. The data are therefore sensitive to optical transient event 



- 7- 



rates A^opt ^ 1 degree ^year ^. 

Unfortunately, not every transient source detected by these searches is sufficiently 
well characterized to say whether or not it could be a GRB afterglow. The High Redshift 
Supernova Search has detected some short duration transient sources without obvious 
host galaxies that might be flare stars in the halo or thick disk of the Galaxy, but might 
potentially be GRB afterglows (Schmidt 1997). Estimating the expected number of flare 
star events is possible (see Garibjanian et al 1990) but beyond the scope of this work. 
Observations at a range of time scales would eliminate this possible source of confusion, 
since stellar flares have characteristic durations of minutes to hours (Krautter 1996), while 
the afterglows of GRB 970508 and 970228 had durations of days. 

Comparing the previously estimated coverage of these surveys to the gamma ray burst 
rate {N^ray ~ 0.01 degree^^ year~^), we see that an absence of optical burst afterglows would 
imply 0-yray/^opticai ^ 0.01. This is already slightly above the ratio fi^ ray /^optical ~ 0.004 
predicted in our illustrative maximal beaming model. On the other hand, detection of even 
one optical transient caused by a gamma ray burster would imply that ray /^optical ~ 0.01, 
i.e., that bursts are strongly beamed. In summary, present optical data cover enough sky 
for a preliminary application of our test, but firm conclusions would probably require 
reanalysis of the data with this test in mind. More definitive results would be possible 
with experiments designed at the outset to search for transients associated with gamma ray 
bursters. 

In the radio, the FIRST survey is using the Very Large Array to map the sky at 
20 cm wavelengths. It is achieving limiting sensitivities of 1 mJy (> 5o") in 3 minutes per 
observation over a field ~ 13' in radius (Becker, White, & Helfand 1995). The observed 
radio flux of the 970508 transient has not yet peaked, but seems likely to peak well above 
1 mJy. Thus, the event rate in radio can be constrained to about the same level as the 
optical rate in a comparable amount of telescope time. Because the difference in event rates 
increases as the wavelength baseline increases in relativistic beaming models, a radio search 
could place the most stringent constraints on gamma ray burst beaming. 

5. Discussion 

In addition to its effect on the transient rate at different wavelengths, beaming of 
gamma ray bursts may have a few other observable consequences. One is that we will 
sometimes be near the center fine of the jet and sometimes at its edge. Yi (1994) pointed 
out that this will broaden the apparent luminosity function of the bursts substantially, and 
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explored the implications of this effect for the statistical properties of the gamma ray burst 
population. 

Another consequence is for the light curve of any individual burst. The predicted 
strength of low energy transients is greatly reduced in extreme beaming scenarios. If the 
burst is beamed into an angle Cm > 1/ro, we can expect a qualitative change in the behavior 
of the transient when the bulk Lorentz factor drops to F < 1/Cm- Before this time, the 
burst will obey the predictions of an isotropic model, while afterwards, a correction factor 
~ (r^^)^ must be applied to all the flux predictions for the isotropic case. 

Moreover, the dynamical evolution of a beamed burst can change qualitatively once 
r < 1/Cm- In ^ spherical model, or while T > 1/Cm, the surface area of the expanding blast 
wave scales as r^. However, the burst cjccta and swcpt-up matter expand at the sound 
speed Cs ~ c/y^ in the comoving frame. The transverse size of the ejecta will thus be the 
larger of r(rn and Cstco, where tco is the time since the burst in a frame comoving with 
the ejecta. The second term dominates for T <C 1/Cm, so that the surface area increases 
more rapidly with r and the ejecta decelerate more abruptly. Specifically, we find that the 
bulk Lorentz factor F oc exp(— r/r^) for 1/Cm ^ T ^ 2, where — [MoTqc^ / (Trpc^)]^^^ . 
(Here Mq and Fq are the initial mass and Lorentz factor of the ejecta, and p is the mass 
density of the ambient medium. Note the implicit scaling Mq oc Cm fixed observable 
gamma ray properties.) The observer frame time also evolves exponentially with r 
in this regime, so that most directly observable properties follow power law evolution, 
with a break when F ~ 1/Cm- Wc explore the consequences in a generalization of the 
Paczynski & Rhoads (1993) afterglow model, whose sole spectral feature is a break where 
synchrotron self-absorption becomes important. We concentrate on the regime where the 
ratio / of swept up mass to initial mass lies between I/Fq and Fq. Under these conditions, 
the observed frequency of peak emission v^^rn, Aux density at that frequency F^^^^^, and 
apparent angular size of the afterglow 9 scale as I'^^m ~ ^e^''^' ^u,®,m ~ ^e^^^^, ^-nd 9 ~ 

^ 3/2 1/2 

before the break, and as i^®,m ~ ^© , -^i/,©,m ~ > and 9 ~ afterwards. Numerical 
integrations confirm these relations, though the transition between the two regimes is 
quite gradual for z/® Combining these scalings with the spectral shape (-FV,© oc u^"^ at 

I/O 

^(B < ^is,m and F^^(Q oc z/® at z/® > z^©,m) yields predictions for the light curve at fixed 
observed frequency. The most dramatic feature is in the light curve shape for > z/©,m, 
which changes from F^,^© ~ ^©^^^ to F^,© ~ t^^. These exponents are generally sensitive 
to the assumed electron energy distribution in the blast wave, and so should be taken as 
illustrative of the type of break expected in a beaming scenario. They could be refined by 
considering models where the electron population is fit to observed data. 

Note finally that the scaling exponents are also sensitive to the burster's environment. 
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As an example, the apparent size in the F > 1/Cm regime becomes ^ ~ %~ if we 

take the medium around the burster to have radial density profile p oc . This implies 
that high resolution studies of these sources (e.g., with very long baseline interferometry) 
would have to follow the evolution of the source over a factor 4 in angular size (and ;^ 10 
in t(^) if they are to distinguish between expansion into a uniform density medium {k = 0) 
and a wind environment (A; = 2). 

To summarize, wc propose a new test for beaming of gamma ray bursts. Observational 
constraints on beaming will help eliminate the dominant remaining uncertainty in the total 
event rate and gamma ray luminosity of the bursts. This in turn will help us determine 
which classes of energetic events have the correct frequency to cause the gamma ray bursts. 
The observations required are feasible with present instruments, and existing data may 
already be sufficient to test the most extreme beaming scenarios. 
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by a Kitt Peak Postdoctoral Fellowship. Kitt Peak National Observatory is part of the 
National Optical Astronomy Observatories, operated by the Association of Universities for 
Research in Astronomy. 
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